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Structural transformations at interfaces are of profound fundamental interest as complex examples of phase 
transitions in low-dimensional systems. They are also of practical importance due to their impact on macroscopic 
mechanical, transport and thermal properties of polycrystals 1 2 . To date, the studies of interfacial phase transitions 
have been focused primarily on "complexions" of intergranular thin films in ceramics 3 5 and pre-melting transitions 
in metallic alloys 6 . Little is known about grain boundary (GB) transformations in single-component metals apart 
from the recently found dislocation pairing transition in low-angle GBs composed of discrete dislocations 7 . Despite 
decades of extensive research, no compelling evidence exists for structural transformations in high-angle GBs. 
Experimental observations by high-resolution transmission electron microscopy (HRTEM) are extremely difficurP^ 
and much of the current knowledge about GB structures comes from computer simulations. In this paper we show 
that the critical impediment to observations of GB phase transformations in atomistic modeling has been rooted 
in inadequate simulation methodology. The proposed new methodology allows variations in atomic density inside 
the GB core and reveals multiple GB phases with different atomic structures. Reversible first-order transformations 
between the GB phases are observed by varying the temperature or injecting point defects into the GB region. 
Due to the presence of multiple metastable phases, GBs can absorb significant amounts of point defects created 
inside the material by processes such as irradiation. We propose a novel mechanism of radiation damage healing in 
metals which may guide further improvements in radiation resistance of metallic materials through GB engineering. 

Recent measurements of diffusion in the Cu £5 (310) GB revealed a marked change in the temperature depen- 
dence of the diffusivity at about 800-850 K^EIl suggesting a structural transformation in this boundary. However, 
the concurrent molecular dynamics (MD) simulations 11 did not detect any structural change other than gradual 
accumulation of disorder and eventually GB pre-melting at high temperatures. Further, the type of structural trans- 
formation suggested by the measurements has not been found in any of the previous simulation studies of this 
boundary. 

Most of the MD simulations conducted to date have employed periodic supercells that did not permit variations 
in atomic density in the GB core, which are possible under experimental conditions. This unphysical constraint 
stabilizes one particular GB structure and prohibits transformations to other structures with different GB densities. 
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We propose a more general simulation approach which facilitates variations in GB density allowing the boundary 
to assume the most thermodynamically favorable structure. This approach reveals the existence of multiple GB 
phases characterized by not only different structures but also different densities. This opens the door, for the first 
time, to direct observations of phase transitions between GB phases with different densities in atomistic simulations 
of metallic systems. 

To demonstrate this new approach, we studied two representative high-angle, high-energy £5 (210) and £5 (310) 
symmetrical tilt GBs in several FCC metals with atomic interactions described by embedded-atom method (EAM) 13 
potentials. The choice of the £5 (310) GB was motivated by the recent diffusion measurements described above^^J 
The long-accepted structure of these £5 boundaries is an array of kite-shaped structural units (Fig.fTJ) and d). This 
structure is easily obtained by joining two perfect crystallites along (210) or (310) planes and applying various 
grain translations followed by static relaxation. It has been suggested, however, that in order to sample all possible 
structures of a GB its density must be also varied by adding or removing atoms ^SEl \\r e applied an algorithm 
which finds equilibrium GB structures at K by varying the GB density along with grain translations and static 
relaxation. In Fig.[T]we plot the obtained energies of the £5 GBs in Cu versus the extra number of atoms relative to 
a perfect crystal plane. The minima in this plot correspond to three different equilibrium structures: (i) the normal 
kites b and d, (ii) split kites c and e having extra 2/5 of the (310) plane or 7/15 of the (210) plane, respectively; 
and (iii) filled kites / having extra 6/7 of the (210) plane. To demonstrate the generality of this finding, similar 
calculations were performed with a different EAM Cu potential and for several other FCC metals including Ag, 
Au and Ni (see Supplementary Information for details). In all cases, the structures shown in Fig. [T] represented 
minima of the GB energy as a function of GB density. Remarkably, and contrary to the common belief, none of 
the tested potentials predicts the normal kites to be the lowest-energy structure of the £5 (210) GB! The split kites 
are energetically more favorable for this boundary and are only slightly higher in energy than the normal kites for 
the £5 (310) GB. 

The effect of temperature on the GB structures was studied by MD simulations with various boundary condi- 
tions starting with the normal-kite structure. First, we used the traditional methodology with periodic boundary 
conditions. In the Cu £5 (310) GB, the kites continued to exist up to ~ 0.75 of the melting point T m = 1327 
K. At higher temperatures the GB begins to disorder and eventually pre-melts near T m (Fig. |2^). This behavior 
known from previous workGS cannot explain the transformation suggested by the diffusion experiments^ 11 * 12 ! A 
more complex behavior was observed in the £5 (210) GB. At about 400 K, the kites were found to transform into 
a structure similar to the filled kites (cf. Fig.[TJ) causing a discontinuity in the GB thickness as a function of temper- 
ature (Fig.[2]3). Similar behavior was found with other EAM potentials for Cu, Au, Ag and Ni (see Supplementary 
Figure 2) suggesting that they represent a generic feature of FCC metals. However, these simulations impose a 
constraint on the supply of atoms into the GB core and do not reveal the full picture of possible structural changes. 

For a more complete investigation, we modified the methodology by letting the GB terminate at an open 
surface. The surface serves as a source or sink of atoms which may diffuse in or out of the GB to automatically 
adjust its density. The results are shown in Fig.[3]for an isothermal anneal at 800 K when the GB meets an open 
surface at one end and is constrained by a wall of fixed atoms at the other. In the simulation of the £5 (310) 
GB initiated in the normal kite structure, a new structure nucleates at the surface and grows inside the boundary. 
This structure can be identified as split kites, whereas the untransformed part of the boundary continues to be 
composed of normal kites. The extra atoms (-40% of the (310) plane) required for the transformation are supplied 
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by diffusion from the surface grove, a process which kinetically controls the transformation rate. Eventually the 
new structure penetrates through and reaches the fixed region. When the fixed atoms on the left are made dynamic 
(thus opening another surface), the entire GB reaches the transformed state (Fig. [3^). 

The stability of the split-kite structure was tested by a series of isothermal anneals. It remained stable and 
well-ordered above -800 K and until approximately 1 150 K when it began to disorder and pre-melt. Upon sub- 
sequent cooling from pre-melting temperatures to 1000 K, the ordered split-kite structure reappeared inside the 
GB. This verifies the thermodynamic stability of the split-kite structure and eliminates any suspicion that its for- 
mation was caused by the presence of the surface grove (an additional proof comes from the point-defect induced 
transformations in the absence of surfaces discussed below). We were unable to demonstrate the reversibility of 
the transformation back to normal kites by cooling the GB below 800 K because of the slow (on the MD time 
scale) diffusion rates at low temperatures. However, the zero-temperature results plotted in Fig.[T]establish that the 
normal kite structure has the lowest energy and will be thermodynamically stable at sufficiently low temperatures. 

Fig. [3]i demonstrates that the old and new structures can coexist and are separated by a line defect with an 
atomic-scale cross-section. This allows us to consider the two GB structures as two-dimensional (2D) phases 
separated a one-dimensional (ID) "phase boundary". The phase transformation is first order and is kinetically 
controlled by GB diffusion. To our knowledge, this is the first observation of coexistence of two GB phases 
separated by a ID "phase boundary" in atomistic simulations. It is also evident that this "phase boundary" is 
associated with a GB sterM The exact nature of this ID defect deserves a separate study in the future. 

Similar transformations were found in the £5 (210) GB starting the simulations from the filled-kite structure 
formed at -400 K. During isothermal anneals at temperatures below 1050 K, the GB transforms into its thermody- 
namically stable phase at low temperatures: the split kites. This new phase grows from the surface and eventually 
penetrates through the entire boundary (Fig. When the boundary with split kites was heated up to 1100 K, 
its structure transformed to filled kites. In this case, however, we were able to observe the reverse transforma- 
tion. Upon subsequent cooling down to 1000 K, the GB structure returned to the split kites, demonstrating full 
reversibility of the phase transformation. Further details of this transformation are shown in Supplementary Figure 
4. At temperatures close to T m this GB pre-melts by disordering of the filled-kite structure as observed earlier^. 
It is interesting to note that the normal kite structure traditionally attributed to this boundary never appears as its 
thermodynamically stable phase. 

During the phase transformations, the GB absorbs/ejects large amounts of atoms from/to the surface. These 
absorption/ejection processes may have important implications for properties of polycrystals under extreme condi- 
tions. For example, in metals subject to radiation by energetic particles, large numbers of vacancies and interstitials 
are formed 21 . Interstitials are much more mobile than vacancies and reach sinks first, whereas the remaining va- 
cancies can form clusters, stacking-fault tetrahedra and other defect complexes degrading the material properties. 
GBs can act as effective sinks of radiation-induced defects and can play a role in the damage healing. Some of 
the mechanisms of defect-GB interactions have been studied by MD simulationsl 22 l 23 l Our work suggests that the 
absorption of point defects can strongly modify the GB structure, causing structural transformations which, in turn, 
can greatly increase the absorption capacity of the boundary. 

The proposed effect is illustrated in Fig. [5] The simulation was intentionally performed in a periodic simulation 
block to exclude all sources/sinks of atoms other than the radiation defects. The GB tested is Cu £5 (310) and its 
initial structure consists of normal kites. This structure remains unchanged during anneals when no atoms are added 
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to or removed from the system. We then introduce 80 interstitials into the GB region, which is equivalent to 40% 
of a (310) atomic plane. This fraction of the plane matches the local minimum of the GB energy corresponding 
to the split-kite structure (cf. Fig. [TJj). Upon anneal at 800 K, the interstitials diffuse into the GB core and get 
absorbed in it initially creating a disordered structure. As the anneal continues, the GB eventually orders into the 
split-kite phase. To demonstrate the reversibility of this radiation-induced transition, we then randomly remove 
the same number of atoms from the GB region, which simulates the insertion of radiation-produced vacancies. 
After a period of disorder the GB structure transforms to back to the normal kites. Thus, in the end of this phase 
transformation cycle the boundary returns to its initial state having annihilated a large amount of radiation-induced 
defects. A similar transformation cycle induced by absorption of point defects was observed in the £5 (210) GB 
(Supplementary Figure 5). 

In summary, we found by atomistic simulations that high-angle metallic GBs can have multiple stable or 
metastable phases with different atomic structures and first-order transitions between them. For the particular 
symmetrical tilt £5 GBs studied here, we found three different phases. The normal kite structure predicted to be the 
K ground state of the £5 (310) GB is consistent with room-temperature HRTEM observations^. The remaining 
two structures have not been seen by HRTEM yet but could have manifested themselves through sudden changes 
in temperature dependencies of GB properties such as diffusivityl im2 l In fact, our simulations show that the GB 
phases found here are characterized by different GB diffusion coefficients and the slopes of their Arrhenius plots 
correlate with the experiment! lm2 l A brief summary of the diffusion simulations can be found in Supplementary 
Figure 6 with more details to be published later. 

All three GB phases found in this work were reproduced in several FCC metals modeled with different atomistic 
potentials, suggesting that they represent a generic feature of FCC crystals. More general GBs may feature a 
larger variety of phases. The existence of GB phases in metals is likely to be a common phenomenon, which 
was previously overlooked due to the overly constrained computer simulation methodology and limitations of 
experimental atomic-level observations of GB structures at elevated temperatures. 

The finding of first-order transitions in elemental high-angle GBs is important as it provides a new testing 
ground for the fundamental theories of GB thermodynamics and GB phase transitions^. Future investigations of 
the exact transformation temperatures, transformation enthalpies as well as the excess thermodynamics properties 
of the ID "phase boundaries" would be of significant fundamental interest. In the future this study could be ex- 
tended to multicomponent systems which may additionally exhibit segregation-induced structural transformations 
with rather complex phase diagrams. In addition to the fundamental questions raised by this work, the present 
simulations suggest a possible novel mechanism of radiation damage healing in which a GB absorbs large amounts 
of point defects by transformation between different phases. The same mechanism of defect absorption can oper- 
ate under other extreme conditions involving high non-equilibrium concentrations of point defects, such as severe 
plastic deformation, diffusional creep, and rapid grain growth in nano-crystalline materials. 
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Figure 1: Results of GB structure and energy calculations at K using an EAM potential for Cu^S. (a) GB energy 
versus atomic density (fraction of atomic plane) for the £5 (310) (left) and £5 (210) (right) GBs. The horizontal 
line marks the energy of the kite structure. The densities corresponding to local energy minima and the respective 
GB structures are indicated. Panels b-f show the GB structures as viewed parallel to the [001] tilt axis (left row) 
and normal to it (right row), (b) and (d) Normal kites, (c) and (e) Split kites, (f) Filled kites. 
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Figure 2: Effect of temperature on the structure of (a) £5 (310) and (b) £5 (210) Cu GBs from MD simulations 
with periodic boundary conditions. The GB thickness was evaluated from the energy density profile across the GB. 
The £5 (310) GB gradually accumulates disorder and pre-melts near the melting point T m , whereas the £5 (210) 
GB undergoes a transformation at 400 K before pre-melting near T m . 
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Figure 3: (a) MD simulation of the Cu £5 (310) GB at 800 K with an open surface allowing variations in GB 
density. The boundary undergoes a first order phase transformation with the split-kite structure nucleating and 
growing from the surface. The normal -kite (b) and split-kite (c) structures are separated by a ID "phase boundary" 
accompanied by a step (d). (e) Completely transformed GB structure in the presence of two open surfaces. 
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Figure 4: MD simulation of the Cu £5 (210) GB at 800 K with an open surface allowing variations in GB density. 
The boundary undergoes a first order phase transformation from filled kites to split kites nucleating and growing 
from the surface. 
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Figure 5: GB phase transformations induced by interstitials and vacancies in a simulation block with periodic 
boundary conditions at T = 800 K. After 80 interstitials were introduced into a 10 A thick layer containing the 
GB, it transforms from the initial normal-kite structure (a) to a disordered state (b) and then to split kites (c). After 
the subsequent introduction of 80 vacancies into the same GB layer, the split-kite structure disorders (d) and then 
transforms back to normal kites (e). The GB states (a) and (e) are identical confirming that the transformation is 
fully reversible. 
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Supplemetary Information 



I GB energies and structures at K 

1.1 Methodology 

The two symmetrical tilt GBs, £5 (310) [001] and £5 (210) [001], were created by joining two crystals rotated 
relative to each other around the [001] direction by the angle of 36.87° and 53.13°, respectively. For each boundary, 
its energy was minimized for several different cross-sections of the periodic simulation block. A cross-section is 
defined by n x m, where n and m are the numbers of lattice periods in the directions parallel and normal to the tilt 
axis, respectively. For the £5 (210) GB we used the following cross-sections: 4 x 2, 6 x 1, 2 x 1, 7 x 2, 9 x 2, 

II x 2, 13 x 2 and 15 x 2. For the £5 (310) GB the cross-sections were 10 x 2, 8 x 2, 9 x 6, 9 x 2, 7 x 2 and 
13 x 2. 

For each GB with a given cross-section, a number of trial configurations were generated and then minimized as 
follows, with the lowest minimized energy reported. First, relative translations of the two grains by vectors within 
a primitive cell of the displacement shift complete (DSC) lattice and all possible, non-equivalent, placements of the 
boundary plane in the direction normal to the boundary were generated. (These different placements are equivalent 
to relative translations by full DSC vectors.) For example, for the 7x2 searches, 3x3x3 vectors within the DSC 
primitive cell were used, and the total number of configurations generated at this stage, including the different 
placements of the boundary plane, was 225. These will be referred to as base configurations below. For some 
searches more points within the DSC primitive cell were used. 

For each base configuration, all atoms within 15 A of the nominal boundary plane were moved 0.001 A in a 
(pseudo) random direction. Atoms that were very near each other were then removed to generate the configurations 
to be minimized. At each step of the removal process, the closest two atoms were both removed and an atom put 
back at the average position of the two that were removed. The process continued until no pair of atoms was within 
about 1/3 of the nearest neighbor distance. This was the first trial configuration to be minimized. After that, one 
atom was removed at a time (from the non-minimized configuration) to generate a configuration to be minimized. 
This was done until no pair of atoms was within about 0.99 of the nearest neighbor distance. In total, about 16,500 
configurations were minimized in a 7 x 2 search for the £5 (310) GB and about 24,000 for the £5 (210) GB. The 
total number N of atoms removed could be larger than the number Nq of atoms in a perfect plane parallel to the 
GB. In such cases, the GB density was determined by dividing N by Nq and expressing the remainder as a fraction 
ofN . 

1.2 Results 

Figure [6] shows several examples of searches with different cross-sections for the £5 (210) GB. For each structure 
corresponding to a local minimum of energy we computed the excess GB energy [U]n, excess volume [V]n and 
two components of the GB stress, rfj and r^ 18 l 26 l The results of the calculations are summarized in Tables W 
and [2] Note that the periodic unit cells of the low energy GB structures are different from each other and can be 
significantly larger than the repeat unit of the lattice. 

The GB structure calculations were also performed for several FCC metals modeled with EAM potentials 
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£5(310) 


Cross-section 


Fraction of (310) plane 


[U] N , J/m 2 


[V] N , A 


r N lt J/m 2 


^,J/m 2 


Normal kites (0 K) 


1 x 1 





0.9047 


0.316 


1.305 


1.774 


Split kites (0 K) 


10 x 2 


2/5 


0.911 


0.233 


0.493 


0.0465 


Split kites (MD) 


18 x 6* 


0.37 


0.920 


0.245 


0.517 


0.0971 



Table 1: Characteristics of different structures of the £5 (310) GB in Cu modeled with the EAM1 25 potential. 
*The MD unit cell represents the size of the simulation block for which the GB properties were calculated and is 
larger than the smallest periodic unit. 



£5(210) 


Cross-section 


Fraction of (210) plane 


[U]n, J/m 2 


[V] N ,A 


t&, J/m 2 


^,J/m 2 


Normal kites (0 K) 


1 x 1 





0.951 


0.322 


1.174 


1.491 


Split kites (0 K) 


15 x 2 


7/15 


0.936 


0.172 


0.297 


-2.12 


Split kites (MD) 


18 x 38 * 


0.46 


0.98 


0.23 


0.446 


-1.276 


Filled kites (0 K) 


7x2 


6/7 


0.953 


0.301 


1.042 


2.299 



Table 2: Characteristics of different structures of the £5 (210) GB in Cu modeled with the EAMlEl potential. 
*The MD unit cell represents the size of the simulation block for which the GB properties were calculated and is 
larger than the smallest periodic unit. 



independently developed by different groups over the past 25 years: EAMl-Cu 25 , EAM2-Cul22l Ag^S, Au^and 
Ni 30 . The results of the calculations for the 7x2 search are summarized in Fig. [7^ and Fig. [7J5 for the £5 (310) 
and £5 (210) GBs, respectively. Note that all the potentials and materials tested exhibit low energy structures 
other than the normal kites at the same GB density (fraction of plane). This demonstrates that the low-energy GB 
structures identified in this work are not specific to the particular Cu potential but are likely to represent a generic 
feature of several FCC metals. 



2 Molecular Dynamic Simulations 

2.1 Methodology 

Molecular Dynamic simulations were performed in the N VT canonical ensemble with a Nose-Hoover thermostat)^ 
using the ITAP molecular dynamics (IMD) program^]. The [001] tilt axis of the GB was aligned parallel to the x 
direction with the GB normal parallel to the y direction (Fig. IS}. The sizes of the simulation blocks used in this 
study are summarized in Table [3] The blocks with periodic boundary conditions (PBC) parallel to the GB where 
used to compute the GB thickness, observe the GB premelting and simulate the GB transformations induced by 
vacancies and interstitials. These blocks had relatively small cross-sections. In other simulations, the periodic 
boundary condition was kept in the x direction, whereas in the z direction the GB terminated at an open surface on 
one side and at a layer of fixed atoms on the other side. In some simulations, the GB terminated at open surfaces 
on both sides. The blocks containing the surface(s) were larger than the periodic blocks. 
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Figure 6: Examples of energy searches for the £5 (210) GB in Cu modeled with the EAM1 potential from Refill 
with different GB cross-sections. The GB energy is measured in the units of 0044, where a is the lattice constant 
and C44 is the elastic constant. Note that the individual searches exhibit sharp minima at different fractions of the 
(210) plane. 
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Figure 7: GB energy calculations for different FCC metals modeled with different EAM potentials: EAMl-Cu^, 
EAM2-Cu^, Ag^, AuE^and Ni^. To facilitate comparison, the GB energy is measured in the units of 0044, where 
a is the lattice constant and C44 is the elastic constant, (a) £5 (310) GB. (b) £5 (210) GB. All calculations used the 
same 7x2 cross-section measured in lattice repeat distances of the bulk crystal. Note that all of these potentials 
identify the normal kite, split kite and filled kite configurations as local minima of the GB energy. However, the 
relative energies of these GB structures are different for different metals. 
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GB 


x, nm 


y, nm 


z, nm 


^ atoms 


(310), pbc 


6.6 


13.0 


6.98 


47,304 


(310), surf 


6.6 


13.0 


35.0 


236,520 


(210), pbc 


6.6 


13.0 


65.8 


46,080 


(210), surf 


6.6 


13.0 


36.2 


253,440 


(210), surf 


9.9 


20.0 


9.6 


143,520 



Table 3: Sizes of the simulation blocks used in this study. The GB was either subject to periodic boundary 
conditions parallel to its plane (pbc) or terminated at one or two surfaces (surf). 
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Figure 8: (a) Simulation block with periodic boundary conditions (pbc) in the x and z directions. In the y direction 
the block terminates at free surfaces, (b) Initial structure of the simulation block with two fixed regions in the y 
direction and one the z direction. One side of the GB in the z direction is terminated at the open surface, (c) GB 
transformation during an isothermal anneal, (d) The GB transformation is completed by creating a second surface 
in the z direction. 



In the blocks with a smaller GB area diffusion-controlled transformations happened faster. Such blocks were 
used for exploratory runs. In most of the simulations with open surfaces, larger block dimensions in the z direction 
were used to ensure that the GB contained large regions inside the bicrystal that were not affected by the surface. 
The GB structures were visualized and analyzed using the visualization software AtomEye developed by Ju 



2.2 MD simulations of GB transformations 

Prior to the MD simulations with open surfaces, the simulation block was annealed at a given temperature with 
periodic boundary conditions parallel to the GB plane. At each temperature, the block was pre-expanded according 
to the thermal expansion factor of bulk copper computed in Ref.^. The grains were terminated with open surfaces 
in the y direction normal to the GB plane to ensure that the normal stress was zero (Fig. [8^). After the anneal with 
PBC, fixed regions were created in the y direction and on one side of the z direction (Fig. [HJ?). The other side of 
the GB in the z direction was terminated at an open surface. At this stage the PBC were maintained only in the x 
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direction parallel to the tilt axis. 

During the subsequent isothermal anneal a GB transformation occurred. The new GB structure nucleated at the 
surface/GB triple junction and propagated inside the bicrystal until the new GB structure reached the fixed region 
(Fig. [8};). At that point the simulation was interrupted and the fixed atoms on the left z side were freed creating a 
second open surface. This resulted in the completion of the transformation (Fig. |SJl). 

The isothermal simulations for the GB transformation were performed at the temperatures of 800 K, 900 
K, 950 K and 1000 K. The MD simulations were also performed with Au and Ag potentials and the expected 
transformations were also observed similarly to the case of EAM Cu 25 . The GB transformation typically occurred 
within the first 5 to 10 ns by nucleation at the surface/GB triple junction. To observe the complete transformation, 
simulation runs as long as 150 ns were implemented. The transformation time depends on the system size as well 
as temperature. 

After the GB transformation has occurred, we cut out an inner portion of the simulation block and recreated 
the PBC in both directions parallel to the GB plane. This procedure gave us a periodic simulation block with a 
uniform GB structure. By comparing the number of atoms in the block with the transformed GB structure and in 
a block with the normal kite structure of the same cross-section, we were able to calculate the change in the GB 
density resulting from the transformation. The obtained GB density is expressed as a fraction of atoms in the (310) 
or (210) plane and is reported in Tables [T] and [2] The numbers obtained correlate well with the results of the K 
calculations. 

The periodic blocks with different GB structures were scaled down to the K lattice constant and relaxed. The 
relaxed structure was used to compute the excess GB energy [U]n, the volume [V]n and the two components of 
the GB stress and 

2.3 Reversibility of the phase transformation in the £5 (210) GB 

The split-kite structure is the most stable phase of the £5 (210) GB in the temperature interval from K to -1050 
K. During an isothermal anneal at 1 100 K the split kites transformed into filled kites. The new GB phase nucleates 
at the surface and propagates into the boundary until the transformation is complete. During this transformation the 
filled-kite and split-kite structures are separated by a step, analogously to the case of the £5 (310) GB discussed 
in the paper. This character of the transformation confirms that it is first order. At the end of the 1 100 K anneal we 
obtained a uniform flat boundary with the filled-kite structure. 

The filled-kite structure was then cooled from 1100 K back to 1000 K and the reverse transformation was 
observed. An the end of the 1000 K anneal the split-kite structure was completely restored. These phase transfor- 
mations are illustrated in Figure [9] 

2.4 GB transformations induced by interstitials and vacancies. 

The GB transformations induced by point defects were modeled in a simulation block with PBC parallel to the GB 
plane. For the £5 (310) GB, the block had the dimensions 6.6 x 13.0 x 6.98 nm 3 and contained 47,304 atoms. 
For the £5 (210) GB, the block had the dimensions 6.6 x 13.0 x 6.58 nm 3 and contained 46,080 atoms. Each GB 
initially had the normal kite structure which was first annealed at 800 K for several ns. Then, 80 and 132 interstitials 
were randomly introduced inside a 1 nm thick layer containing the £5 (310) and £5 (210) GBs, respectively. These 
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Figure 9: High temperature transformations in the £5 (210) GB terminated by two open surfaces. The initial 
split-kite structure taken from a prior anneal at 1000 K transforms into filled kites at 1100 K. Upon subsequent 
cooling and anneal at 1000 K, it returns back to the split kites. Intermediate stages show the two GB structures: 
one growing from the surface and the other disappearing. It is concluded that the phase transition temperature is 
between 1000 K and 1100 K. 
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numbers of atoms correspond to -40% and -50% of the (310) and (210) planes, respectively. After the atoms were 
inserted, the isothermal anneal was continued at 800 K. The initial kite (respectively, filled kite) structures quickly 
disappeared and significant atomic rearrangements took place in the boundary layer. However, within 1 to 100 ps 
of the simulation the boundaries re-ordered into the split-kite structure (Fig.[T0|). 

The obtained split-kite structures of both GBs were then taken as initial states to introduce vacancies in order to 
demonstrate transformations into metastable phases. The vacancies were introduced by randomly removing atoms 
within a 1 nm thick layer containing the GB. The number of vacancies was equal to the number of interstitials 
introduced previously. After 1 to 100 ps of the isothermal anneal at 800 K and a period of disorder, we observed a 
transformation into the normal kite structure in the £5 (310) GB and the filled-kite structure in the £5 (210) GB. 
Figure ( 10 1 demonstrates these transformations for the £5 (210) GB. The corresponding images for the £5 (310) 
GB appear in the paper. 



2.5 GB diffusion in the £5 (310) GB 

As discussed in the main text, experimental measurements of GB diffusion in Cu bicrystals containing the £5 (310) 
GB have been made by Budke et al.^and more recently by Divinsky et al.^l The measurements revealed a marked 
change in the slope of the Arrhenius plot of the diffusivity at about 800-850 Kj 1 1 * 12 l To test the prediction that this 
change is caused by a structural phase transformation in this boundary, we performed MD simulations of the 
GB diffusivity using the methodology proposed in Ref r 5 " 37 l The simulations were performed separately for the 
normal-kite and split-kite structures. Recall that the normal kites are stable at low temperatures and transform to 
split kites at temperatures 800 K and higher. 

The results of the diffusion calculations are shown on the Arrhenius diagram in Fig. [TT] Note that the two 
GB structures have different diffusivities as well as activation energies. Furthermore, it is observed that the high- 
temperature split-kite phase has a lower activation energy (smaller slope) than the low-temperature normal-kite 
phase (larger slope). This trend closely correlates with the experimental measurements ■LU121 as d oes the approxi- 
mate temperature at which the change in the slope of the Arrhenius diagram is expected. This agreement strongly 
suggests that the change in the GB diffusivity measured in the experiments^^] was caused by the phase transfor- 
mation in this boundary. Further details of the diffusion simulations will be published elsewhere. 
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Figure 10: Phase transformations in the £5 (210) GB induced by interstitials and vacancies in a simulation block 
with PBC at T = 800 K. The GB transforms from filled kites (a) to split kites (c) after 132 interstitials were 
introduced in a 1 nm thick layer containing the GB. The GB transforms from the split kites (c) back to filled kites (e) 
after 132 vacancies were introduced in same layer, (b) and (d) demonstrate that significant atomic rearrangements 
take place during the transformations. Note that the initial (a) and final (e) structures are identical. This metastable 
filled-kite structure does not spontaneously transform into the more stable split-kite structure due to the absence of 
sources and sinks of atoms in the PBC block. 
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Figure 11: Arrhenius diagram of the diffusion coefficient D calculated for two structures of the £5 (310) GB in 
Cu. Note the difference in the activation energy between the two GB structures. 
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